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ABSTRACT 

We have used the Optical Redshift Survey (ORS; Santiago et al. 19 
due to fluctuations in the galaxy density field out to distances of 8000 km s^^. At large scales where 
linear theory applies, the comparison of this gravity field with the observed peculiar velocity field offers 
a powerful cosmological probe, because the predicted flow field is proportional to the parameter ri°-^/6, 
where VL is the matter density and b is the bias of the galaxy distribution. The more densely sampled 
ORS gravity field, t o excellent approximation, matches that of the earlier IRAS 1.2-Jy redshift survey 
(Fisher et al. 1995a), provided (3 is reduced by a factor bopt/bjRAS ~ 1-4. Apart from this scaling, the 
most signihcant din'erence between the ORS and IRAS fields is induced by differing estimates of the 
over-density of the Virgo cluster. Neith er of these gravity f ields is consistent with the peculiar velocity 
field constructed from the full Mark HI (Willick et al. 1997a) sample. We find that a simple but plausible 
non-linear bias algorithm for the galaxy distribution relative to the mass has a negligible effect on the 
derived fields. We conclude that the substitution of optical for IRAS catalogues cannot alone resolve 
the discrepancies between the IRAS gravity field and the Mark HI peculiar velocity field. 

Subject headings: cosmology — dark matter — galaxies: clustering — large-scale structure of universe 



1. INTRODUCTION 

With the advent of large, uniform redshift surveys of 
galaxies, increasingly stringent tests of models for the 
large-scale str ucture of t he u niverse have become pq ssi- 
ble (see, e.g., |Dekel 1994| and ^trauss fc WiUick 1995| for 
reviews). Under the assumption that large-scale flows of 
galaxies in the universe are a response to the underly- 
ing distribution of matter, peculiar velocity measurements 
are a critical probe of cosmology and large-scale structure. 
Redshift surveys allow one to construct an estimate of the 
gravity field based on the observed distribution of galax- 
ies. Using linear or quasi-linear gravitational instability 
theory, one can predict the velocity field from this gravity 
field, a prediction which is approximately proportional to 
the parameter f3 = f{il)/b, where b is the galaxy bias, f 
is a f unction which is well approximated by ( Peebles 
1980 ), and H. is the matter density. 

Until recently, the Infrared Astronomical Satellite 
(IRAS) 1.2-Jy flux-limited redshift survey was the only 
available nearly full-sky catalogue. The original redshift 
survey of IRAS-sclected galaxies was limited to a sam- 
ple of 2658 galaxies w ith 60-micron flux greater 1.936 Jy 
(Strauss et al. 1992a), but this was later extended to a 



sample of 5321 galaxies brighter than 1.2 Jy (Fisher et al 



1995a). A sample of 13,000 galaxies compl ete to 0.6 Jy is 



nearing completion (Canavezes et al. 199g) and is eagerly 



awaited. The principal advantage of these JRAS-selected 
samples is their insensitivity to extinction within the Milky 
Way, allowing extremely large sky coverage, 87.6% (11.01 
sr) for the 1.2-Jy sample, with a single, linear instrument. 
Nearly full-sky coverage is essential for estimation of the 
peculiar gravity field. 

However, it is well known that the distribution of IRAS 
galaxies is biased with respect t o the overall distribution 
of galaxies ( Strauss et al. 1992b|) , primarily because IRAS 
under-counts the dust-free early-type galaxies which con- 
gregate in cluster centers. The IRAS surveys are also more 
dilute than is desirable, typically including only 1/3 of the 
known spiral galaxies in volumes where the selection func- 
tion is high. Optically-selected surveys can eliminate this 
bias, and deeper redshift surveys will overcome the shot 
noise problem associated with the sparse sampling of the 
IRAS survey. The recently compl eted Optical Redshift 
Survey (ORS; Santiago et al. 1995 ) is currently the best 
available catalogue for our peculiar velocity analysis. The 
ORS is a concatenation of three optically-selected sam- 
ples covering most of the sky with |6| > 20°, and it is the 
best current approximation to a full-sky optically-selected 
catalogue. Until a redshift catalogue is constructed from 
the 2MASS survey ([Stiening, Skrutski, fc Capps 1995| ), the 
ORS is likely to remain the closest competitor to the IRAS 
full-sky catalogue. Previous work (Freudhng, da Costa, & 
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THE ORS-PREDICTED VELOCITY FIELD 



Pellegrini 1994) suggests that the IRAS and optical grav- 
ity fields are consistent, but with the superior data now 
available, it is clearly of interest to determine the nature 
of any deviations between the gravity fields of the optical 
and IRAS samples, as these could substantially affect com- 
parisons to observe d peculiar velocity fields. Direct density 
field comparisons ( Santiago fc Strauss 1992 ; Strauss et alj 
1992b) between the ORS and other redshift su rveys will 
be carried out in detail by Strauss et al. (1998) . 

Comparisons of the gravity field derived from the IRAS 
survey with observed peculiar velocity samples have been 
carri ed out using a number of complementary methods 
(see Strauss fc Willick 199| for a review), with more re- 



cent s tudies perforra e d by Davis, Nusser, fc Willick (1996 



herea fter, DNW96)|, IWillick et al. (1997b)|, |Ricss et al 



The Optical Redshift Survey (ORS) is a combination 
of three individual declination-limited catalogues: the Up- 
psala General Catalogue (UGC) at northern declinations 
{S > —2.5°), the European Southern Observatory (ESO) 
catalogue in the south {S < —17.5°), and the Extension 
to the Southern Observatory Catalogue (ESGC) in the re- 
maining strip just south of the Celestial Equator. The 
sample selection a nd galaxy distributio n have been de- 
scribed in detail (Santiago et al. 1995), and the selec- 
tion functions for the galaxies in each sam ple have been 
carefully calibrated (Santiago et al. 1996). We use the 
flux-limited versions of the UGC and ESO catalogues; for 
the ESGC catalogue, only a diameter-limited sample was 
available. All velocities have been converted to the Lo- 



1997M, ISigad et al (1998)|, |da Costa et al. (1998)|, and 



Willick fc Strauss (1998) These studies show that the 



cal Group (LG) frame using the transformation of Yahil 
Tammann fc Sandage (1977)| . 



JffAS-predicted flow field is at least qualitatively very sim- 
ilar to the measured velocity field, lending strong support 
to the gravitational instability model for the growth of 
large-scale structure. Although the qualitative alignment 
of the JffAS-predicted and Mark III peculiar velocity fields 
is remarkably good, detailed quantitative comparisons are, 
at present, less satisfactory. DNW96 find large coherent 
residuals between the IRAS and Mark III fields, primarily 
in the dipole component at lar ge scales, which preclu de a 
conclusive determination of f3. Willick et al. (1997b) find 
better alignment of the IRAS gravity and Mark III velocity 
fields using the rather different VELMOD technique, but 
this analysis is limited to redshifts within 3000 km s^^- 



Hudson (1993a, 1993b) constructed a full-sky map of 
the local density field from optically-selected galaxies, us- 
ing the UGC and ESO catalogues limited to a redshift of 



8000 km 



Complete redshift information did not exist 



Willick fc Strauss (1998) find good agreement on larger 



scales, but conclude that this requires a re-calibration of 
the TuUy-Fisher relations for the various subsamples going 
into the Mar k HI catalogue incon sistent with the published 
calibr ation (Willick et al. 1997a ). On the other hand, ^ 
Costa et al. (1998) find a better match between the large- 
scale IRAS gravity field and the SFI I-band sa mple of spi 



ral fie ld gal axy peculiar velocity measurementsj Giovanell: 
' 19971). The SFI sam ple uses the same Mathewson 



et al. 

Ford 

Mark III, but the transformation of the data to a com- 
mon system differs from the Mark III treatment. A recent 



fc Buchhorn (1992)| data for the Southern sky as in 



for these catalogues at the time of his analysis, and his sky 
coverage was limited to 67%. Furthermore, his procedure 
for merging the UGC and ESO into a single catalogue may 
have led to spurious asymmetries between the North and 
South. The complete ORS redshift sample allows us to 
overcome many of these limitations. 

Our analysis relies on a spherical harmonic expansion 
of the density field and therefore requires data of uniform 
quality over the entire sky. Unfortunately, there are several 
limitations which prevent the ORS from being a uniform, 
all-sky survey. Because the ORS galaxies are optically 
selected, the catalogue is affected much more adversely 
by extinction than is the IRAS survey. The ORS does 
not cover the strip at Galactic latitudes |6| < 20°, while 
the IRAS survey excludes only the strip with |6| < 5°, 
two small unobserved strips at higher latitudes, and small 
confusion-limited regions (about 12.4% of the sky in total). 
In addition, we only use ORS data in regions where the 
Burstein fc Heiles (1984)| blue extinction is low: Ab < 0.7 



summary of the current confusion is given by Davis (1998) 



Might some of the discrepancy disappear if we were to use 
an optical redshift survey to construct an estimate of the 
gravity field? 

Thi s paper addresses the question of the construction 



magnitudes. Finally, there are 64 plates, each 5° on a 
sid e, which are missing from the ESO-LV survey (Fig. 5a 
of Santiago et al. 1995| ). The accepted ORS data cover 
about 62% of the sky, and the completeness of the cata- 
logues within these regions is excellent. 

To produce an all-sky catalogue, we use the IRAS sur- 
vey to fill in the areas in which ORS data are absent or 



the ex tinction is high. The IRAS catalogue (Fisher et al 



of a three-dimensional gravity field based on the ORS and 1995a) contains 5321 galaxies covering 87.6% of the sky 



IRAS redshift catalogues. We construct our estimate of 
the gravit y field using the redsh ift-space procedure de- 
scribed by Nusser fc Davis (1994). Section 2 describes the 
galaxy catalogues and our procedure for merging them. 
The method for obtaining the gravity and velocity fields 
is briefly outlined in §3, with details of the spherical har- 
monic formalism in the Appendix; more details have been 
presented elsewhere. Section 4 describes the differences 
between the ORS field and the field derived from IRAS 
alone, including a brief discussion of alternatives to the 
simple linear biasing scheme, and §5 discusses the main 
conclusions. 



2. THE GALAXY CATALOGUES 



We use a version of this catalogue which Yahil et al. (1991) 



have carefully augmented with fake galaxies in the Zone 
of Avoidance (and other excluded regions) without biasing 
the statistics of the galaxy distribution. This all-sky IRAS 
catalogue contains 6010 galaxies. The IRAS data are only 
used in regions where there are no accepted ORS data, so 
each of the four surveys cover disjoint regions on the sky. 

The ORS survey has the advantage that it is much 
more densely sampled than the IRAS survey within cz ~ 



8000 km 



This allows the local gravity and velocity 



fields to be constructed with much smaller sampling er- 
rors. However, the ORS selection function drops exponen- 
tially with distance, while the IRAS galaxy distribution 
has a long power-law tail extending to higher redshifts. 
Beyond 8000 km s~^, the ORS data become unacceptably 
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sparse. We therefore use only IRAS data over the whole 
sky for redshifts cz > 8000 km s^^, extending the analysis 
to cz = 18,000 km s"^ 

To generate a density field in redshift space, we first 
merge the galaxy catalogues into a single, all-sky, uniform 
catalogue. Each catalogue j has a selection function 
such that for homogeneously distributed galaxies, the ex- 
pected galaxy density at a point s = {s,9,ip) is nj(f>j{s), 
where fij is the true average density of galaxies in the cat- 
alo gue volume. Since w e are given the selection functions 
0j ( Santiago et al. 1996| ) , we calculate the density in (some 
part of) the volume of catalogue j using 



by 



1 



1 



(1) 



where the sum is over galaxies i in the volume of inter- 
est. The sky coverag e Vtj of each of the surveys is known 
(Santiago et al. 1995), so 



V- — -O s'^ 



(2) 



is known. Galaxies with > Smax,j are not included 
in the analysis. We calculate the density within Smax — 
8000 km s^i for the three ORS catalogues and the IRAS 
catalogue, although, as mentioned above, the IRAS analy- 
sis is carried out to 18, 000 km s~^. The densities we com- 
pute are for objects sufficiently luminous to be included in 
a volume-limited sample out to radius cz — 500 km s~^. 

The fractional density fluctuation is given by the ra- 
tio of the observed galaxy counts to the number expected 
given the selection function. This can be calculated by 
assigning each galaxy i a weight [nj(j)j{si)]~^ . In our 
merged catalogue, we multiply these weights by factors 
which normalize the density in each catalogue to the den- 
sity of IRAS galaxies within the accepted volume of the 
optical subsample. This allows each ORS subsample to 
have a density which, due to inhomogeneities, is not nec- 
essarily the same as the overall density, but it matches 
the different ORS subsamples to the uniformly-selected 
IRAS catalogue. Thus we calculate ujhasj, the density 
of IRAS galaxies within the volume of catalogue j, rela- 
tive to niBAS, the overall density of IRAS galaxies. The 
weight for any galaxy is then given by 



1 



niRAS./ 



nj0j(sj) niRAS 



(3) 



This weighting ensures that the density of the merged cat- 
alogue remains approximately constant across catalogue 
boundaries. S ome parameters of th e catalogues are listed 
in Table (cf. [Santiago et al. 19961 ). 



3. CALCULATION OF THE VELOCITY FIELD 

Nusser & Davis (1994) show that in linear theory, the 
peculiar velocity field in redshift space is irrotational and 
can therefore be derived from a potential: v(s) = — V$(s). 
If the angular dependences of the potential and galaxy den- 
sity field = 6p^/p^ (both measured in redshift space) 
are expanded in terms of spherical harmonics as described 
in Appendix A, the potential is related to the density field 



ds 



Im 



ds 



1 i{i + 1)$/™ 



l + (3 



_ld\n 



(3 s 



s din s ds 



(4) 



This is a modified Poisson equation, where Sf^ is the spher- 
ical harmonic coefficient of the redshift-space galaxy den- 
sity field, s is the redshift, and (3 was defined above. The 
last term on the right hand side corrects for the "rocket 
effect" , which arises because we weight each galaxy by the 
selection function (j) evaluated at its redshift rather than 
its (unknown) distance. 

We now briefiy summarize our computationally efficient 
method for solving the above equation. We first compute 
a smoothed density field on a spherical grid in redshift 
space. We construct a 32 x 32 angular grid equally spaced 
in Galactic longitude and latitude, with 52 bins in redshift 
out to 18, 000 km s^^, providing more than adequate res- 
olution given the smoothing of the fields. The separation 
of the redshift bins increases in proportion to the mean 
IRAS inter-particle spacing, {n(j))~^^^, in order to reflect 
the decreased sensitivity of the surveys at high redshift. 
The Gaussian-smoothed galaxy density field at grid point 
n is given by 



l + <5ns«) 



1 



(27r)3/2a3 



^(s„-s,) /2< 



(5) 



where the sum is over all galaxies in the merged catalogue. 
The Gaussian smoothing width for the cell, (T„, is equal 
to the mean IRAS inter-particle spacing at that redshift 
(or 100 km s~^ when the inter-particle spacing is smaller 
than this). The increased smoothing at larger redshifts is 
essential to prevent divergence of the increasing shot noise; 
our procedure is qu alitatively close to the optimal Wiene r 
filtering procedure (Lahav et al. 1994; Fisher et al. 1995"b| ), 
and the result is a signal-to-noise ratio in the density field 
that is roughly constant over radial bins. A correction 
factor must also be applied because part of the Gaussian- 
weighted volume centered on a cell will fall outside the 
survey cutoff radius. 

The spherical harmonic coefficients of the density field 
are calculated on radial (redshift) shells. Once these co- 
efficients have been computed for the merged cata- 
logue, the modified Poisson equation for the potential is 
solved. The velocity field can then easily be constructed 
from V = — V$. 

Our analysis is done purely in redshift space, assuming 
linear theory and a one-to-one mapping between distance 
and redshift along any given line of sight. These assump- 
tions become invalid in clusters, where highly nonlinear, 
virializcd motions ca use a radial distort ions into "fingers 
of God." We follow [Yahil ct al. (1991)| in collapsing the 
galaxies in the six nearest clusters of their Table 2 to a 
single redshift. 

A limitation of the density field constructed by merg- 
ing disparate catalogues is that we can only use one se- 
lection function for the dln^/ds term on the right hand 
side of the modified Poisson equation (see Appendix B 
for details). The maps below were constructed using the 
IRAS selection function, which, however, is not a very 
good approximation to the ORS selection functions for 
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IRAS 


38% 


6010 


3020 


0.0488 0.0488 


1.047 


UGC 


34% 


3246 


2992 


0.0799 0.0473 


1.015 


ESO 


19% 


2412 


2187 


0.109 0.0361 


0.775 


ESGC 


9% 


1351 


1203 


0.195 0.0570 


1.223 



Table 1 

Numbers of galaxies, densities, and catalogue weights for the catalogues used in the ORS density field analysis. The 
IRAS values include only the portion of the survey outside the accepted ORS volumes; the IRAS sky coverage is 100% 
for cz > 8000 km s^^. The full-sky IRAS density for cz < 8000 km s^^ (optimal for minimum-variance estimates of the 
density) is uihas = 0.0466 Mpc"^. 



cz > 8000 km s~^. To gauge the importance of the differ- 
ing catalogue selection functions in the rocket effect correc- 
tion, we have compared the ORS peculiar flow field using 
both (j^uGC and (jiiRAS ■ (Differences among the three ORS 
selection functions are negligible in comparison to the dif- 
ference between ORS and IRAS.) Although the rocket ef- 
fect correction is substantial over the redshifts of interest, 
the 4>uGC and 4>iras flow fields are identical (to within 
^ 30 km s~^, with negligible systematic differences over 
the sky) for cz < 8000 km s~^, where the selection func- 
tions behave similarly. At greater distances, however, the 
fields diverge dramatically, particularly in the dipole. Be- 
cause we use only IRAS data for cz > 8000 km s^^, we 
are fortunately able to avoid this complication. 

4. COMPARISON OF THE FIELDS 

A detailed statistical comparison of the ORS and IRAS 
gravity fields is difficult due to the limited sky coverage of 
the ORS. Because we require an all-sky catalogue for the 
spherical harmonic analysis, we have filled in the Galac- 
tic plane with IRAS galaxies; the two fields are therefore 
not independent over the whole sky. In addition, the ORS 
shares some galaxies in common with the IRAS survey. 
Nevertheless, the qualitative comparison of the two fields 
can tell us whether or not improved rcdshift surveys might 
provide a better match to the Mark III data than DNW96 
found for the IRAS survey. Detailed Wiener-filtered re- 
constructions of the IRAS den sity and gravity field in real 
space have b een carried out ( Webster, Lahav & Fishei 
1997 ; see also Strauss & Willick 1995); our analysis, how- 
ever, is done entirely in rcdshift space. 

4.1. ORS and IRAS Density Fields 

The density field contours on several radial shells are 
illustrated in Figure ^, where it is apparent that the ORS 
and IRAS fields are qua litatively very simil ar (compare 
with the similar figures in [Santiago et al. 1995] ). The JRAS 
survey's under-counting of cluster centers is most evident 
on the near (cz = 1000 km s~^) radial slice, which con- 
tains Virgo {e = 284°, 6 = 75°), Ursa Major {£ = 145°, 
b = 66°), and Fornax {£ = 237°, b = -54°). With the 
Gaussian smoothing for this shell {a — 345 km s~^), the 
ORS over-density for Virgo is 6^ ~ 10, compared with 
5^ « 6.7 in the IRAS field (recall that the smoothing in- 
creases with radius, so density contrasts appear smaller on 
more distant shells). Note that if the IRAS galaxies dif- 
fered from the ORS galaxies only by an extra bias factor 
of 2/3, which would match the peak Virgo over-densities 



for the two fields, then we would expect that a compari- 
son of observed peculiar velocities with the ORS gravity 
field would yield a value of (3 lower than the corresponding 
number for the IRAS gravity field by the same factor. 

On the next redshift slice {cz — 3000 km s^^), the 
Hydra-Centaurus complex is clearly visible around = 
300°, b = 20°). The maximum over-densities are ap- 
proximately 25% smaller for IRAS than for ORS. How- 
ever, the smoothing on this shell is somewhat larger 
(fj — 550 km s^^), and the clusters are quite near the ORS 
b = 20° cutoff. The linear structure of the Supergalactic 
Plane is visible near £ = 135°, and there are pronounced 
voids (Jf w -0.9) towards £ = 265°, b = -50° and I = 5°, 
5= 10°. 

At cz — 5000 km s^^, the Perseus-Pisces region around 
{I — 140°, b = —25°) is clearly much better sampled by 
the ORS survey than by the IRAS survey, but the quali- 
tative features of the two fields are identical. The smooth- 
ing at this redshift is 815 km s~^. The b < —20° parts 
of the Pavo-Indus-Telescopium region are also visible near 
I — 330° . The density contrasts on more distant shells are 
much less pronounced, and beyond 8000 km s~^ the two 
fields are identical by construction. 

4.2. Predicted Flow Fields 

Figure || shows the IRAS- and ORS-predicted radial pe- 
culiar velocity fields (Local Group frame) on redshift shells 
for j3 = 0.6 (note these are not the same as the redshift 
shells in Fig. H). Again, the overall agreement of the fields 
is quite good. The nearby shell (cz = 500 km s~^) is dom- 
inated by Virgo-centric infall, with galaxies towards Virgo 
and galaxies in the opposite direction of the sky flowing 
away from us. The ORS-predicted flow towards Virgo 
is approximately twice as large as the IRAS-predicted 
flow for the same value of (3. There are also substantial 
{v ^ 300 km s^^) predicted flows towards us from the Su- 
pergalactic poles, but the differences in these regions be- 
tween the ORS and IRAS predictions are more modest 
than for the flow toward Virgo. 

At cz = 2000 km s~^, the dipole component of the flow 
begins to dominate as a result of the motion of the Local 
Group, with infall on the back side of Virgo and outward 
flow in the opposite direction. The flow of galaxies away 
from us as it falls into Hydra-Centaurus is modest but visi- 
ble at this redshift. The dipole pattern grows with redshift, 
as evident in Figure ^c, and except for an overall scaling 
and slight change of orientation, the differences between 
the ORS and IRAS fields diminish. 
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The monopole, dipole, and quadrupole terms as a func- 
tion of radius for the ORS- and JRAS-predicted flow fields 
are shown in Figures for the same value /3 = 0.6, 
the high end of the "most hkely" range of DNW96. The 
quantity plo tted is Qmi as defi ned in Appendix A. As em- 
phasized by [Davis (1998)| and [Nusser fc Davis (1994)1 the 
monopole and dipole terms at a given rcdshift arc sensitive 
only to the interior mass distribution, which should be well 
sampled. Tidal fields introduced by errors in the density 
field of the dilutely sampled exterior region will only affect 
the quadrupole and higher multipoles. Statistically signif- 
icant disagreement between the measured and predicted 
velocity field for the dipole component is thus a sign that 
something is amiss, either a systematic error in one or 
both of the catalogues, incorrect treatment of an impor- 
tant component of non-linear fiow, or non-uniform relative 
bias in the galaxy distribution. These figures compare the 
multipole moments of the IRAS- and ORS-predicted ve- 
locity fields, to see whether the DNW96 discrepancy con- 
tinues to hold for the ORS data. 

The ORS-predicted dipole terms (Fig. P are all signif- 
icantly larger than the IRAS-predicted terms. The differ- 
ences in the Galactic X (m = 1) and Y {m = —1) com- 
ponents arise primarily behind Perseus-Pisces, but the Z 
{ni = 0) difference is mainly due to Virgo. The quadrupole 
terms agree quite well, except that the m = — 1 magnitude 
is somewhat lower in IRAS, and the ORS data resolve the 
m = peak around Perseus-Pisces a bit more sharply. We 
should point out that because each of the ORS subsam- 
ples are normalized to the IRAS density field (eq. ||^), the 
IRAS and ORS dipoles are not as statistically independent 
as we would otherwise like them to be. 

Figures |^-|| also show the results when the ORS (3 value 
is lowered from 0.6 to 0.4. Note that equation ^ is only 
linear in /? in the limit P 0. For P of order unity, red- 
shift distortions lead to more complex behavior; thus the 
P = 0.6 and P — 0.4 ORS multipole terms are not simply 
proportional to each other. The overall amplitude of the 
ORS-predicted dipole for P = 0.4 agrees extremely well 
with the JRAS-prcdicted dipole for P = 0.6, but the in- 
creased weight of Virgo enhances the relative strength of 
the Z-component. Lowering the value of P has a fairly 
modest effect on the quadrupole terms, although the peak 
in the overall quadrupole amplitude at 2500 km is sig- 
nificantly reduced. The m = —1 component is brought 
into excellent agreement with the IRAS prediction, and 
the m = ±2 amplitudes are lowered somewhat below the 
IRAS prediction. In any case, the incorporation of the 
ORS data has worked in the wrong direction to correct 
the mAS-Mark III dipole residual of DNW96, because 
the Y component of the Mark III field is much too small 
relative to the IRAS prediction, and the ORS data only 
drive these components further out of agreement. 

The dipole of the IRAS-predicted velocity field does not 
seem to converge to the CMB dipol e (£ = 268°, b = 27° 



in the LG frame) particularly well (Strauss et al. 1992c; 
Webster et al. 1997); it consistently points about 30° away. 
With their larger Virgo-centric infall, the ORS data only 
exacerbate this situation. At cz = 6000 km s~^, the ORS 
dipole points to ^ = 250°, b — 55°, about 15° higher in lat- 
itude than the IRAS dipole. This is reflected in the larger 
magnitude of the I — 1, m = multipole component of the 
field (Fig. H), which points in the Galactic Z direction. 



4.3. Predictions for Mark III Galaxies 

The predicted peculiar velocities for the Mark Illgalax- 
ies in three rcdshift slices are plotted in Figures The 
IRAS plots differ from those in DNW96 only in the projec- 
tion used and the orientation of the longitude coordinate. 
Note that the IRAS predictions use /3 — 0.6 while ORS 
uses P — 0.4, consistent with the apparent relative bias 
of the two fields o f a factor of 1.5 (Strauss et al. 1992b; 
Hermit ct al. 1996). In general, the difference between the 
ORS- and JRAS-predicted fields is quite small. In fact, 
the magnitude of these residuals is generally lower than 
the residuals seen by DNW96 in their tests on mock cat- 
alogues. The strong velocity shear seen in the Mark III 
sample across the Hydra-Centaurus complex is plainly ab- 
sent from either of the rcdshift surveys; indeed, this region 
shows only very small differences between the fields. The 
largest systematic differences are seen for Perseus-Pisces; 
our peculiar motion away from this complex is smaller in 
the ORS-predicted field. The detailed statistical analysis 
of DNW96 shows discrepancies between IRAS and Mark 
III which would not be alleviated by our ORS gravity field. 

4.4. Effects of Non-linear Bias 

In general, the bias of the galaxy distribution is the 
function which relates fiuctuations in the galaxy number 
density Sp^ to fluctuations in the overall matter distribu- 
tion Sp. The functional form depends on the details of 
galaxy formation, which are poorly understood. For sim- 
plicity, it is usually taken to be a linear function indepen- 
dent of scale, so that the bias param eter b has a single 
numerical value. Semi-analytic work (Kauffman, Nusser 



& Stcinmetz 1997; see also Weinberg 1995) suggests that 
most physical bias mechanisms lead to a nearly linear bias 
on large scales, but the value of the bias depends quite 
strongly on the power spectrum, galaxy type, and galaxy 
luminosity in a given sample. This is of course problematic 
for flux-limited surveys, in which the most distant galaxies 
are also the most luminous. The ORS data in fact show 



some evide nce of a weakly scale-dependent bias (Hermit 
et al. 1996| ). 



Provided that the bias is linear, our technique will re- 
cover the value of /3 « il'^-^/b. Non-linear bias in the 
galaxy distribution will introduce systematic errors in the 
analysis. We investigate the impact of a non-linear bias 
by postulating a gridded mass density fluctuation field 
5 — ^B{d^) — yS] /b, where B is some non-linear function. 
The mean value B over the set of grid points must be 
subtracted because the fluctuation field has zero mean by 
definition. 

Most prescriptions for non-linear bias are functions of 
the local over-density S. With the variable smoothing of 
our computed density field, however, it is difficult to imple- 
ment such a non-linear bias in a spatially uniform fashion. 
But the mean value of the field is J = independent of 
smoothing, and it is therefore possible to construct a non- 
linear bias using 6 — as a pivot. A simple non- linear 
form is the piecewise linear function 



BiSS) 



{2/3)63 



(6) 



as shown in Figure ^. This in effect makes galaxies more 
strongly clustered than dark matter in over-dense regions, 
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and less strongly clustered in under-dense regions, relative 
to the case of linear bias. In other words, we are stealing 
dark matter from the clusters, filling the voids with more 
than would be expected on the basis of the galaxy distribu- 
tion. The galaxy voids are thus substantially deeper than 
the matter voids for reasonable values of b (Cen & Ostriker 
1992, 1993). This is plausible (though admittedly ad hoc) 
because galaxy formation is expected to be more efficient 
in dense reg ions, and it is consistent with the results of 



simulations (White et al. 1987; Dekel 1997; Lemson et al_ 
19980 We emphasize that our crude piecewise linear bias 
is only meant to be illustrative; a more complex algorithm 
could not readily be implemented because of the variable 
smoothing of our density field. 

Even the fairly dramatic "knee" we have introduced in 
the bias at = has quite a modest effect on the pre- 
dicted flow field. The overall amplitude of the flow is re- 
duced, as expected because the dark matter distribution 
is now less clustered. However, the magnitude of this re- 
duction is only 10% (see Fig. ||), and the best-fit value 
of /3 would therefore only increase slightly. The direction 
of the dipole is altered by only ~ 2°, and other qualitative 
features of the flow field are unchanged. It seems quite 
unlikely that even a rather extreme uniform bias prescrip- 
tion could plausibly eliminate the coherent residuals seen 
between Mark III and our predicted velocity fields. 

5. CONCLUSIONS 

The gravity field derived from the ORS is quite similar 
to the field derived from the IRAS survey, provided the 
value of f3 is multiplied by a factor of « 2/3, roughly con- 
sistent with the relative strength of clustering in the two 
surveys. Some differences from this simple scaling are ob- 
served; most significantly, the dipole of the ORS-predicted 
flow field is rotated towards Virgo and is further than the 
IRAS prediction from the CMB dipole. This non-scaling 
behavior arises because the ratio of ORS to IRAS den- 
sity fluctuations is larger in the local Supercluster than in 
other regions. 

DNW96 have found systematic discrepancies between 
the Mark III and JfJAS-predicted flow fields, particularly 
a large coherent dipole residual, which cannot be physical 
because the dipole field depends only on the interior mass 
distribution. Estimating the gravity field from the ORS 
survey has provided a check on the influence of the known 
shortcomings of the IRAS catalogue (namely, its sparse 
sampling and under-representation of cluster cores) . Since 
the IRAS and ORS gravity flelds are so similar, we con- 



clude that the redshift-space catalogues are not themselves 
the source of the DNW96 discrepancies. 

Compared to IRAS, the ORS gravity field will gener- 
ally push us towards lower values of /? in order to fit a 
given sample of peculiar velocity measurements. The re- 
sult of DNW96 for the JRAS-Mark III comparison would 
then suggest /3 <^ 0.45, but this comparison with Mark 
HI remains suspect due to to the substantial quantita- 
tive discrepancies betw een the IRAS and Mark III flelds. 



da Costa et al. (1998)| , on the other hand, showed that 
the SFI sample agrees quite well with the IRAS gravity 
fleld for P = 0.6; presumably this sample would also be 
consistent with the ORS field for (3 « 0.4. Another seri- 
ous concern for the Mark III comparison is the fact that 
complementary methods for attacking the velocity field 
problem have not converged to a con sistent, unambigu- 
ous result. The POTEN T analysis of ^igad et al. (1998)| 
and Dekel et al. (1993)| favors larger values (/3 ^ 1) than 
the ITF and VELMOD analyses. These methods weight 
the data in very different ways, and given that the fields 
show some level of inconsistency, it is not surprising that 
different analyses of the same data yield different answers. 

The sources of the DNW96 discrepancies between the 
IRAS gravity field and the Mark HI peculiar velocity mea- 
surements remain uncertain. These discrepancies might be 
telling us that galaxies cluster in a way which is not closely 
related to the clustering of the mass, or that the linear 
gravitational instability model provides an inadequate de- 
scription of flows on the scales of interest. Alternatively, 
there may be unknown system atic biases in the Mark HI 
data; see also the discussion in Willick & Strauss (1998) 



We have investigated a simple but plausible generalization 
of the usual linear bias model for galaxy clustering, but 
this does not signiflcantly affect features in the predicted 
velocity fleld; even a rather sizeable change in the slope of 
the bias relation at (5^ = causes only a modest reduction 
in the amplitude of t he fleld. Imp roved full-sky peculiar 
velocity surveys (e.g., Strauss 1997) of larger galaxy sam- 
ples are eagerly awaited. Such surveys are quite a difficult 
undertaking, but they may go some way towards clarifying 
our understanding of large-scale flows in the universe. 
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APPENDIX 



SPHERICAL HARMONIC FORMALISM 

We solve for the linear peculiar velocity fleld from the density fleld expanded in spherical harmonics on radial (redshift) 
shells. The conventional spherical harmonics arc deflned to be 



where 



'21 + 1 {I -my: 



An {l + my. 



1/2 



(Al) 



(A2) 



and the P/m's are the associated Legendre polynomials. Since the density and velocity flelds are real, it is conveni ent to 
work with the real-valued spherical harmonic functions, which differ slightly from their complex cousins. Following Bunn 
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(1996), we define them as 

( 1 TO = 

yirn = {-Y'KiraPlrn{cOSe) x\ \/2 COS TO(^ TO > (A3) 

[ V2 sin|TO|(^ TO < 0. 
The coefficients aim of an arbitrary real field A(s) are related to the field by 

aim{s) = JdnyimA{s) (A4) 

OQ I 

A{s) = '^'Myi'n- (A5) 

It is convenient to normalize the coefficients of multipole order I by the value of yio at the North Galactic Pole; we 
therefore define 

Qm = ( 4^ J (^^) 

With this normalization, cn is the amplitude of the dipole vector projected along the Galactic X axis {£ ^ b = 0°), ci^_i 
is the amplitude along the Y axis {£ = 90°, b = 0°), and cio is the amplitude along the Z axis (6 = 90°). The amplitude 
of the dipole is the quadrature sum of the three coefficients. The quadrupole terms of the real-valued spherical harmonics 
are 

3^20 = iy|(3cos2 0-l) (A7) 
1 / 15 

3^21 = -W — sin 26* cos (A8) 

4 V TT 

1 / 15 

3^2-1 = -a/ — sin 26* sin (A9) 

4 V TT 
1 /T5 

3^22 = -\/ — sin2 0cos2(/5 (AlO) 

4 V TT 

3^2.-2 = 7W — sin2 0sin2(^, (All) 

4 V TT 

where Galactic longitude is ^ = <^ and latitude is 5 = 90° — 6. 

ROCKET EFFECT CORRECTIONS 

A limitation of constructing a merged density field out of different catalogues and then sol ving for the v elocity field 
is that equation (^ only allows for a single radial selection function in the rocket effect term (Kaiser 1987). This term 
corrects for the fact that, when constructing the density field, we have evaluated the selection function at each g alaxy's 
redshift rath er than its (unknown) distance. This can be a substantial correction to the redshift-space density field (Fishei 
et al. 1995b| ). 

For the ORS survey, unfortunately, the selection function </)(s) is not independent of angle on the sky, but depends on 



the cat alogue and the extinction along the line of sight. The linear equation for the velocity potential $ of Nusser & Davis 
is 



(1994) 



where / « f2°-^, b is the bias, and is the galaxy over-density in redshift space. A general selection function which 
is a function of position on the sky [9, ip) will introduce couplings between different modes of the spherical harmonic 
expansion. Expanding <I> on each redshift shell as 

oc / 

'^lra{s)Yira{e,^) (B2) 

i=0 m=-l 



and similarly for (5^, multiplying through by and integrating over solid angle, we find that the term subtracted from 
5^ in equation (jj) becomes 

Ul'rn' [ d(cOs6>) dip y^m^'m' a!"'^ ' (B3) 



^ 00 / 



i'=0 m' = -l' 



s fr^ , ,, ./(47r) olns 
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where u = ~d^/ds is the radial pecuhar velocity. For a selection function independent of azimuthal angle ip, this 
expression reduces to 




This expression replaces the — (u/„j/s)(o?ln0/c?lns) term in equation (Q). For our catalogues, the angular dependence of 
the selection function is approximately just a step discontinuity at |6| = 20°, so by symmetry, even (odd) values of I' 
make no contribution to odd (even) values of I. The utility of the spherical harmonic expansion breaks down unless the 
I' = I term is dominant. For the analysis at cz < 8000 km s~^, the corrections are negligible because the ORS and IRAS 
selection functions track each other quite well. A deeper analysis, or one in which the selection functions of the merged 
catalogues are not so similar, may need to take account of these corrections. 
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Fig. 1. — Smoothed density field contours on radial shells and catalogue galaxies in redshift slices. The Gaussian 
smoothing width subtends approximately 10° on each shell. The long-dashed line shows the ORS |5| = 20° cutoff. All 
redshifts are measured in the Local Group frame, (a) The nearby shell cz = 1000 km s^^ with galaxies in the range 
< cz < 2000 km s^^, contour spacing 1.0, Gaussian smoothing width a = 345 km s~^. (b) The intermediate shell 
cz = 3000 km s~^ with galaxies in the range 2000 < cz < 4000 km s~^, contour spacing 0.5, Gaussian smoothing width 
(7 = 550 km s~^. (c) The distant shell 5000 km s~^ with galaxies in the range 4000 <cz< 6000 km s~^, contour spacing 
0.5, Gaussian smoothing width a = 815 km s~^. 
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Fig. 2.— Radial slices of the ORS- and IRAS-predicted peculiar velocity fields in the Local Group frame for /? = 0.6. 
Dot-dashed contours show negative (approaching the Local Group) velocities, (a) cz — 500 km s^^, contour spacing 
50 km s~^. (b) cz = 2000 km s~^, contour spacing 200 km s~^. (c) cz = 5000 km s~^, contour spacing 200 km s~^. 
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Fig. 3. — Monopole coefficient cqo of tlie real-valued spherical harmonic expansion of the ORS- and JRAS-predicted 
peculiar velocity fields. The solid line is the JRAS-predicted field for /3 = 0.6, the dashed line is the ORS-predicted field 
for the same /?, and the dotted line is the ORS-predicted field when /3 is reduced to 0.4. 
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Fig. 4. — As in Figure H, but for the dipole coefficients cim- The quadrature sum is shown in the lower right panel. 




Fig. 5. — As in Figured, but for the quadrupole coefficients C2m- The quadrature sum is shown in the lower right panel. 
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Fig. 6. — Predicted peculiar velocities for the Mark III galaxies based on ORS (/3 = 0.4, top) and IRAS {(3 — 0.6, middle) 
for the near redshift slice {cz < 2000 km s^^). Circles indicate positive peculiar velocities (directed away from the LG), 
and crosses indicate negative velocities. The bottom plot shows the difference between the ORS and IRAS predictions. 




Fig. 7.— As in Figure ^, for the middle redshift slice (2000 < cz < 4000 km s ^). 




Fig. 8.— As in Figure ^, for the distant redshift slice (4000 < cz < 6000 km s^^). 
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Fig. 9. — (a) A simple non-linear form for the bias function B (solid line), consistent with our variable smoothing 
procedure. The fractional fluctuation in the galaxy density is 6^. (b) Effect of the non-linear bias above on the dipole 
amplitude of the ORS peculiar flow field with /3 = 0.6. The dashed line is for linear bias. 



